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Abstract. We present a microscopic description of edge excitations in the quan- 
tum Hall effect which is analogous to Feynman's theory of superfluids. Analytic 
expressions for the excitation energies are derived in finite dots. Our predictions 
are in excellent agreement with the results of a recent numerical diagonalization. 
In the large N limit the dispersion law is proportional to qlog^. For short range 
interactions the energy instead behaves as . The same results are also derived 
using hydrodynamic theory of incompressible liquids. 
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In the last few years a considerable interest has been devoted to the study 
of the edge excitations of 2D charged systems (quantum dots) in strong magnetic 
fields [1-5]. These studies are motivated by the fact that edge excitations, differ- 
ently from bulk modes, are gapless and are consequently particularly relevant for 
the thermodynamic behavior at low temperature. The macroscopic picture under- 
lying many of the available theoretical works is the hydrodynamic description of an 
incompressible 2D liquid characterized by the propagation of edge waves with drift 
velocity v = where E is the electric field generated by the electrons at the edge 
of the droplet (see, for example, ref.[4]). According to this picture the dispersion 
should be 

M E 

where M is the angular momentum carried by the wave and R oc y/N is the radius 
of the dot {N is the number of electrons). The above picture presents however a 
serious difficulty since the electric field generated by 2D charged clusters exhibits 
a logaritmic enhancement at the border. Moreover considerable theoretical effort 
has been recently devoted to models where electrons interact through effective short 
range forces [5] and where consequently the concept of electric field cannot be used. 

The purpose of this work is to give an answer to this problem and to provide 
an explicit formula for the dispersion law employing either a microscopic approach 
based on Feynman's theory and a macroscopic description based on classical hy- 
drodynamics. We will also discuss the crucial role played by the neutralizing back- 
ground in ensuring the stability of the system. 

In the following we will consider the Hamiltonian 
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where the vector potential has the form = \yB^Ay = —^xB, the second term 
is the e-e Coulomb interaction, while the last term accounts for the Coulomb in- 
teraction with the neutralizing background for which we make the simple choice 
Pion{f) = Po0(^ ~ R)- Due to charge neutrality the electron density in the interior 
of the dot coincides with the ion density po- The radius R is then fixed by the 
normalization condition pqttR^ — N. 

We will focus on the case of integer filling (u = 1) where it is possible to derive 
important results in an analytic way. This fixes a relationship between the density 
Po and the magnetic field. We also assume that the magnetic field is strong enough 
that mixing of states in higher Landau levels by the Coulomb interaction can be 
neglected. An important consequence of the integer filling is that the ground state 
of the Hamiltonian (2) is a Slater determinant [6] built up with the single particle 
states 

where I = 0,1, N — 1 and £ = {hc/eBy/"^ is the magnetic length. This Slater 
determinant carries angular momentum Lq = N{N — l)/2 and is the state with the 
lowest energy because of the presence of the confining potential. In the absence of 
such a field the electrons would in fact occupy states with higher values of I in order 
to minimize the Coulomb repulsion. 

The electron density corresponding to the ground state is characterized, for 
large N, by a constant value po = interior of a circle of radius R = £y/2N 

and by an edge thickness of the order of the magnetic length £. The electric field 
generated by the electrons can be easily calculated for large N. At the border of 
the cluster we find the result E{R) = \og{f3N) with (3 = 16.4. This equation 
explicitly shows the anticipated logaritmic enhancement. 
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Our approach to the study of the dispersion of the edge excitations is based 
on the idea, currently considered in the hterature [5-7], that the lowest edge state 
I M > carrying angular momentum L — Lq+M is naturally excited by the collective 
operator 

N 



si = 2-"/^ Y.Ci, - = E f^4.MC. (4) 

k 1 

where q and c| are the electron annihilation and creation operators relative to the 
single particle states (3). The operator (4) corresponds to the projection of the 
usual multipole operator onto the lowest Landau level. It is the edge analog of the 
projected density operator used in ref.[8] to study bulk excitations in the fractional 
quantum Hall effect. 

The key point of the work consists of the ansatz | M >= | >. This ansatz 
has the form of the approximation employed by Feynman to describe the density 
excitations of superfluid '^He. The same method has been succesfuUy applied to 
study bulk excitations in the fractional quantum Hall effect [8]. 

For large N the transition density ptriT^) =< M \ 5{r — r^) | > associated 
with the state | M > takes the typical form ptr{^) p'{r)exp{—i9M) of an edge 
wave with wave vector q = M/ R. 

The excitation energy of the "Feynman" state can be written in the following 
form: 

< Q|[^M,[g,^l,]]|Q> 

<0| [Sm,sI] |0> 



< M 


H 


M > 


< M 


M > 



<M) = - <o\H\o >= ^ (5) 



where we have used the fact that the operator Sm annihilates the ground state. For 
a fixed value of M this bound is expected to coincide with the exact dispersion law 
in the large N limit. 
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In the following we calculate the numerator of eq.(5) using the Hamiltonian 
(2). The kinetic term does not contribute to the the commutator since S\j does 
not excite higher Landau levels. In order to calculate the double commutator we 
find it convenient to use the Tamm-Dancoff relation (see, for example, ref.[9]) 

<o\[Sm,[h,sI]]\o>= Y.s*^,a (6) 

minj 

where Armnj = Sijdmniem - e^) + and = Y^fjo^ ^mjmj + ^mm are the 

single particle Hartree-Fock energies written in terms of the usual two-body and 
one-body matrix elements V^^^^ and V^^. The quantity Smi — 8M,m-i{N^Y^'^ is 
the matrix element of the operator Sm between the ground state and the particle- 
hole state cJ^Q I >. It is worth noticing that result (6) is an exact one since the 
ground state is a Slater determinant. 

Let us first calculate the contribution to the double commutator (6) arising from 
the e-e interaction. By using non trivial relationships involving matrix elements of 
the e-e interaction, we find ioi N » M the most important result 

< I [Sm. S^mW I >= M2[Fa,(M) - Fn{1)] (7) 

whith the quantity Fn{M) defined by 

Fn{M) = [ dridr2'ilJ*N+M{T^iWN-M{^2)-, rV'Jv(ri)V'Jv(r2) . (8) 

J I I"! — r2 I 

Due to the 1/r behavior of the Coulomb force the function Fn{M) exhibits a 




logN dependence. However the logN term vanishes in the difference (7) which 



then behaves as N 



The contribution of the confining potential can be also calculated in the large 
N limit. With our choice for pion 

it is possible to relate this contribution to the 
electric field generated by the electrons and we find 



< I [Sm, [V-*, SU I >= N^M \ /off(/3iV) (10) 



Equations (7-10), together with the result < | [Sm, S]^] \ >= MN^, allow 
us to calculate the excitation energy of the edge modes through the Feynman's 
relation (5). The energy can be written as the sum of the two contributions: 



e = 1= > — r; e""^* = ^=log{BN) . (11) 

7reV2N ^^2l-V 27rey/2N ^ ^ ^ 

The quantity e®~® is always negative and vanishes in the M = 1 mode as a conse- 
quence of the translational invariance of the e-e interaction. For M = 2 it coincides 
with the result recently found in ref.[6]. Notice that the sum of the two terms (11) 
must be positive in order to ensure the stability of the ground state. 

An important property revealed by eq.(ll) is that the e-e term does not de- 
pends linearly on M in contradiction with the eq.(l). Actually for large M we 
find e^^(M) = --M^log{^M) with 7 = 0.96. This logaritmic behavior is the 
consequence of the Coulomb interaction as revealed by eqs.(7-9). Use of the short 
range force aV^5(rj — r^) yields a different dependence: 

- ^^.N-^'m - M^) (12) 

Our predictions for e^~^ and e*^ turn out to be in excellent agreement with the 
results obtained in ref. [5] through a numerical diagonalization of the e-e interaction 
(see figure 1). In fig. 2 one clearly sees the convergency of the results of ref. [5] to 
the asymptotic expressions (11) and (12) holding for N » M. 
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For large M it is natural to introduce the parametrization M = qR = q£v2N 
and the dispersion law, given by the sum of the two contributions (11) takes the 
form 

2 

e{q) = + e^^* = —qlog^ (13) 

TT q 

with go = \/^/\/2^7 = 3.0/£. Note that in the sum (13) the logN terms arising 
from e^~^ and e*^ cancel each other. 

For the short range force the dispersion law instead becomes 

e''^(?) = -:r^5'- (14) 



The applicability of both results (13) and (14) is fixed by the long wavelength 
condition q£ « 1. 

In the last part of the work we show that results (13-14) can be directly ob- 
tained using classical hydrodynamic theory. We treat our system as a charged 
incompressible liquid {diw = 0) characterized by a velocity potential of the form 
X = ;^oe*^^~'*~*'^* where x and y are the directions orthogonal and parallel to the 
border respectively (the liquid occupies the plane x > 0). This implies Vy = —ivx 
and the Euler equation ^ — — -|-vxB — -r— Vp' takes the form 

^ at mc mpo ^ 

p' = —i'mpo{uj + ujc)— ■ (15) 

In eq.(15) Uc = eB/mc is the cyclotron frquency and p' is the oscillating part of the 
2D pressure. The excess of pressure p' is compensated at the border by a restoring 
force produced by the electric field. This force can be calculated by varying the 
energy of the electric field U = I/Stt / E'^dV with respect to the displacement 
^{y) of the border of the liquid. Thus the boundary condition takes the form 
p'^^o ~ ~^U/5^. The calculation of U is drastically simplified if one notices that 
the main contribution to the integral originates from the region q^^ < r < q~^ 
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where r is the distance from the border and Qq is a cutoff length of the order of 
the widht of the border. In this region the deformed edge can be considerd as a wire 
of charge density e^{y)pQ generating the electric field E — 2e^{y)po/r. Integration 
of the electric energy in the xz plane then gives 

U = e^pllog^ J e{y)dy (16) 

showing in a clear way the physical origin of the logaritmic term. Since the x- 
component of the velocity of the fluid at the border is given by = d^/dt, we 
finally obtain the boundary condition p^^^^ = —2ie^ p^^log^ . The Euler equation 
(15) then gives rise to the dispersion law 

oj{oj + 0Jc) = 2—poqlog— . (17) 

2 

Equation (17) yields, for low q and a;, the result e(q) = u^qlog^ where we have 
used the relation po = ^'^^ defining the filling factor u. When u = 1 the hydrody- 
namic formula coincides with the microscopic result (13). This qlog^ dependence 
has been alreday derived in ref. [10-11]. 

It is finally interesting to recover result (14) for short range forces. In this case 
the relevant restoring force originates from the surface energy U — ^ J { ^^gy ^ )^dy 
where a is the surface tension. The resulting dispersion becomes: 

u;(u} + ujc) = -^q^ (18) 
mpo 

yielding, after identifying a = —apo/7r\/27r£^, result (14) in the low q regime. Note 
that in the figures we have reported the results of ref. [5] with a < corresponding 
to negative excitation energies. 

The exact equivalence between the microscopic and macroscopic results dis- 
cussed in this work confirms in a clear way the general statement that electrons 
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in a strong magnetic field exhibit a behavior typical of Bose superfluids and that 
consequently their dynamics is properly decribed by the equations of classical hy- 
drodynamics [12]. Both in the Coulomb and short range cases the dispersion law 
however differs from the linear law e = qv currently considered in the literature and 
reveals new interesting features exhibited by edge excitations in the quantum Hall 
effect. 

We wish to thank R. Ferrari for many useful discussions and R.L. Schult for 
providing us with his numerical results. L.P. likes to thank the hospitality of the 
Department of Physics and the financial support from the Centre ECT* at the 
University of Trento. 
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FIGURE CAPTIONS 

Fig.l. Lowest excitation energies e®~® and e^'^ for N = 400. Crosses are taken 
from ref.[5], while dots are the predictions of eqs.(ll) and (12). The units are e^/l 
and 10~^a/Z'^ for the Coulomb and short-range interactions respectively. 

Fig. 2. A'^-dependence of the lowest excitation energies e*^"*^ and e*'' for M = 10 
Crosses are taken from ref . [5] , while the dashed lines are the asymptotic predictions 
(11) and (12). For the units see fig.l. 
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